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Abstract 

We derive new explicit bounds for the total variation distance between two convolution 
products of n S N probability distributions, one of which having identical convolution factors. 
Approximations by finite signed measures of arbitrary order are considered as well. We are 
interested in bounds with magic factors, i.e. roughly speaking n also appears in the denomina- 
tor. Special emphasis is given to the approximation by the n-fold convolution of the arithmetic 
mean of the distributions under consideration. As an application, we consider the multinomial 
approximation of the generalized mult inom i al dis tribution. It turns out that here the order of 
some bounds given in lRood (|2001r ) and lLohl (|l992l ) can significantly be improved. In particular, 
it follows that a dimension factor can be dropped. Moreover, better accuracy is achieved in the 
context of symmetric distributions with finite support. In the course of proof, we use a basic 
Banach algebra tec hnique fo r meas ures on a measurable Abelian group. Though this method 
was already used bv lLe CamI (|196Cll ). our central arguments seem to be new. We also derive new 
smoothness bounds for convolutions of probability distributions, which might be of independent 
interest. 
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ation distance. 



MSC 2000 Subject Classification: Primary 60F05; secondary 60G50, 62E17. 



1 Introduction 

1.1 Aim of the paper 

Approximations of distributions of sums of independent random variables are needed in nearly all 
branches of probability theory and statistics. Many results for normal and compound Poisson ap- 
proximations are nowadays available. However, if the distributions of the summands are similar to 
each other, much better accuracy can be achieved using identical convolutions of a certain distribu- 
tion. In the present paper, we give total variation bounds for the accuracy of such approximations 
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in a general framework, i.e. for probability distributions on a measurable Abelian group. We also 
consider higher order approximations by finite signed measures. All bounds contain magic factors, 
i.e. roughly speaking n appears in the denominator. As a consequence, this en ables us to d erive 
mult i dimen sional results, some of which improve the order of bounds obtained in iRoosI (|200ll ) and 



Loh 



(|l992l ). It should be mentioned that Loh used Stein's method in a more general situation of 
dependent random variables. However, it seems to be unclear, whether Stein's method can be used 
to reproduce the results of the present paper. Furthermore, it turns out that our bounds have a 
better order in the case of symmetric probability distributions with finite support. Our proofs are 
ba sed on a combi nation of some Banach algebra related techniques, which in principle were used 
bv lLe CamI (| 19601 ). On the other hand, the core arguments given in Sections 14. II and 14.21 seem to be 
new. Further, the smoothness estimates for convolutions of probability distributions in Section [4.11 
might be of independent interest; for instance, see (j35p and (j38p . 

We n ote th a t, at the beginning of our investigation, we tried to improve one of the central 
results of IRoosI (j200ll ). see ([5]) and discussion thereafter. But unfortunately we were not able to 



use the multidimensional expansion of that paper for any substantial improvement. Surprisingly 
it turned out that it is better to forget the dimension, so to speak, and to use the properties of 
measures on a measurable Abelian group. This should explain, why we use this somewhat abstract 
approach. 

The paper is structured as follows: The following two subsections are devoted to the notation 
and a review of known results. In Section [21 we present and discuss our main results. To get a first 
impression of the results of this paper, the reader may consult (fTSj) . (fT6]) . and ([TH]) . In Section [3l 
we give some numerical examples. The proofs are contained in Section [H 



1.2 Notation 

Let (jC, +,A) be a measurable Abelian group, that is, +) is a commutative group with identity 
element and ^ is a c-algebra of subsets of X such that the mapping {x,y) i-^ x — y from 
{X X X, A ® A) to (j£, A) is measurable. We note that it is more convenient to formulate our 
results in terms of distributions or signed measures rather than in terms of random variables. Let 
!F (resp. M.) be the set of all probability distributions (resp. finite signed measures) on {X,A). 
Products and powers of finite signed measures in M are defined in the convolution sense, that is, 
for V,W & M. and ^ G ^, we write ^^^(A) = J^V{A — x) dW{x). Empty products and powers of 
signed measures in M are understood to be I := /qi where Ix is the Dirac measure at point x G X. 
Let V = — V~ denote the Hahn- Jordan decomposition oiV & M and let \V\ = + V~ be its 
total variation measure. The total variation norm of V is defined by ||y|| = |y|(jC). We note that, 
in the literature, often the total variation distance sup^g_4 [^(A) — G(A)| = W\F — G|| between 
F,G £ J- is used. In this paper, however, all distances will be given only in the total variation 
norm. With the usual operations of real scalar multiplication, addition, together with convolution 
and the total variation norm, is a real commutative Banach algebra with unity I. For V £ Ai 
and a power series g{z) = X]m=o^™-^™' i'^^n £ ^) converging absolutely for each complex z £ C 
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with \z\ ^ ll^ll, we define g{V) = Ylm=o (^mV"^- The above assumptions imply that the hmit exists 
and is an element of the Banach algebra Ai. On the other hand, the definition of g{V) can also be 
understood setwise. The exponential of y G is defined by the finite signed measure 

oo ^ 

e^ = exp(y) = e M. 

m=0 

We note that e^ is not necessarily a non-negative measure. Further, exp(t(F — /)) is the compound 
Poisson distribution with parameters t G [0, oo) and F ^ T . \iF and G are non-negative measures 
on (X, ^) and F is absolutely continuous with respect to G, we write F <^ G. For F £ and 
A £ A, F\^ is the restriction of F to the set A. The complement of ^ G ^ is denoted by A^. Set 
= and n = {1, . . . ,n} for n G N = {1,2,. . . }; further, for n G Z+ = NU{0}, set Uq = {0, . . . ,n}. 
For a set J, let \J\ be the number of its elements. For x G M, let [xj = sup{n G Z | n ^ x} and 
[x] = inf{n G Z I n ^ x}. Always, let O'' = 1, 1/0 = oo, and, for k £ Z, Yl'^=k = be the empty 
sum and 11^^ = 1 t^e empty product. For a G C and b G Z+, let (^) = nm=i(^ — m + l)/m. For 
a, 6 G M, set a A 5 = min{a, b}. 



1.3 Known results 



We first discuss some important results for discrete distributions on X = M"^, {d G N) with the usual 
addition. Let 

d d 
Ho = 1, Hr = /e,, (r G d), Fj = ^pj^rHr, U G n, n G N), F = ^PrHr, (1) 

r=0 r=0 

where, for r G dp, pj^r ^ [0, 1] with Ylt=oPj,r = ^iPr = "-"^ J2]=iPj,r > Oj and G M'^, (r / 0) is 
the vector with 1 at positi o n r a nd otherwise. 

y. Theorem 1 and Lemma 2) proved with the help of Stein's method 



In the case d = 1. 



Ehm (1991 



that the total variation distance between the Bernoulli convolution rij=i ^^'^ binomial law F 
can be estimated by 



72 

62 



mm 



1 ^ 
L npiPQ } II 



^ 272 min 



nPiPo J 



(2) 



where 7^ = X]^=i(Pi ~Pj,i)^j ^ N). Her e, the estimates depen d on the behavior of the so-called 

jll992)), and on the closeness of all pj^, 



Barbour et al. 



magic factor {np^po)^ (cf. Introduction in 
(j S n), which is reflected by 72. In Theorem 3 of iRoosI ((20001), a Krawtchouk expansion was used 
to show that an absolute constant C > exists such that, if 72 > 0, then 



^ C6' minU, 



I73l 



+ 



1 



l2\fnpiPQ npipo 



72 



nPlPo 



For example, it easily follows that || 11^=1 — -F" || ~ Y^2/(7re) 9 as —>■ and npipQ — > 00. 
Here, ~ means that the quotie nt of both sides tends t o one. Further results in this and a more 
general context can be found in ICekanavicius and RoosI (120061 ) and the papers cited there. 
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The multivariate case d G N was investigated by iLohl (|l992l ) using Stein's method. He gave 
an estimate for the closeness between the generahzed multinomial distribution rij'=i 
multinomial distribution F". This bound contains certain functions Ci,C2 > of p^, (r G dg), 
which can be estimated from above by absolute constants, if all Pj.'s are uniformly bounded away 
from and 1. In his Theorem 5, he showed that, if n ^ 2 and max{Cin~^/^, C2[2(n — 1)]"^} ^ 1, 
then 

n n 

1 1 JJ - F" ^ 2 ^ \Pj,riPr2 - Pj,r2Pv, I ^ri ,r2 , (3) 

j=l j=10^ri<r2^d 



where 



'ri,r2 



C2 . /2(n-l) 
m ' 



n — 1 



Co 



+ 



Co 



2(n - 1] 



2 2Ci 
+ 



mm 



n 



n n 

{ll^^-m,n),ll(i-m,r2)}. 



1=1 



i=l 



The quantities Ci,C2 can be given explicitly as 



Ci= sup [Ci(ri,r2) A Ci(r2,ri)], €2= sup C2{ri,r2,rs) 

0^ri<r2<d 0^ri,r2,r3^Q!: 



where, for ri,r2,r3 G dg, 
Ci(ri,r2) = 



1 



1/2 



+ 



C2(ri,r2,r3) 



2 3 

Pri Pra ^^Przi'^ - Pr 
(2,2 

' J_ I ^ 2 _ , ^ 2 _ 

2,3, 1 



1/2 



(4) 



+ 



J_ , J_ , 1 _ 

Pri Prg ep^gCl-Prs) 



1/2- 



if r2 = rs, 



if r2 / rs. 



If d = 1, then it follows from Ehm's result and the equality '^o<^ri<r2^d\Pj,riPr2 ~ Pj,r2Pri\ — 
\Pi ~ Pj.il' Loh's bound is not of the best possible order, because of the exponent of \'pi — pj^i\ 
and the logarithmic term. It turned gut that a b ound better than ^ can be given using a 
multivariate Krawtchouk expansion, see iRoosI (|200ll . Theorem 2, Corollary 1). Indeed, 



r=l 



where 



2-v^ 



^ 10.15 and 



A sometimes more precise bound is 



4 1 

Kf) = y^{Pr -Pj,r)^min{-, h {r £ d). 

^ e npj,pQ ) 



1 - E?.i 



r=l 



(5) 



(6) 



In contrast to Q, for d = 1, the bounds in ([S]) and ([5]) have the same order. In the general case, 
from dS]) and Cauchy's inequality, it follows that || YYj=i ~ ^|| ^ C^dY^f^^ We note that 
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this estimate is of the same accuracy as ([5]) when the S{r), (r G d) are of similar magnitude. In 
view of this bound, one might wonder, whether the dimension factor d can be dropped. However, 
as shown in iRood (|200ll . Remark 2 after Proposition 2), this is not generally possible. But if we 
concentrate on the estimate with the magic factors, i.e. 

j=l r=l j=l •^'^•^0 

the more general results of this paper imply that C3 d can indeed be replaced by the constant 21.88, 
see Example 12.11 below. It should be mentioned that here the Hy, {r E dg) need not just be the 
Dirac measures as in ([T]). 



2 Main results 

In what follows, we present bounds which are small when the Fj E (j E n) are close or when n 
is large and the Fj are not too different. Our first result is the following. 

Theorem 2.1 Zei n E N, Fi, . . . , F„, G E .F, Fq = F = i E"=i Fj, 



JCn:\J\=kj&J k=0 



For j,k e n and m E N, set M,-,fc = {Fj - G)GL("^^)/*=J , 
Set 



Y^'^=i \\Mj^k\\^, and 



rie,a = max 

k e n\e 



1 / ? 



+ 1^ 



k,2 



, Eng, a E [0,00)), r]^ = ripQ, 



Ci = sup 

xe(o,oo) 



ln(2 - (1 - xy 



0.694025 . . . 



(a) Let a E [0,oo), i E Uq, and [3 = \a{l + l)/2] . //r?^ < (2eCl)"^ then 



(8) 



In particular, for a = 0, we have 



(2eci7?,)(^+i)/^ 
1 - v^2e Ci ■ 



(9) 



(b) Assume that, for each j E ng, Bj E .4 exists suc/i that Fj|^c ^ G and let fj denote 



Radon- Nikodym density of Fj\^c with respect to G. For i E tiq' then have 



1/2 



■E[^7|^(l|(^.-G)|Bj| + |(F,-G)(i?,)l)+Vr( / (/.-l)'dG)"" ".(10) 



1/2 
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We note that, if G = F, then = and r]^ ^ simplifies to rj^^ = max^jg^y might ask 

why we gave the complicated estimate ([8]). However, it turns out that, in special situations, the 
order of rj^ ^ for a > can be much better than that of rj^. See Proposition 12.11 below involving a 
bound for ry^ i instead of just the estimate (jlOp . Further, the reason why we formulated Theorem 12. II 
in its present general form without the assumption that G = F is given with Lemma 14.31 and 
Example 14.21 below. 

Let us first discuss the simple case when a = 0. 

Remark 2.1 Let the assumptions of Theorem 12.11 hold. In what follows, whenever we consider Vk 
or Wk for a specified number k G Z^, we assume that k ^ n. 

(a) For ken, let Ffc = y,'J_. (G-F,)K We have Vq = I,Vi = n(F-G), V2 = ^{n^(F-G)^ -T2), 



and, similarly as in iRoosI (|200d . formula (10)), it can be shown that 

fc-i 

k 



^ k-l 

Vk = —rY,^j'^k-j, (ken). 



j=0 

This formula can easily be used to evaluate the signed measures Wi for a given i. In particular, 
we have Wq = G" and 

Wi = G'' + n(F-G)G'"\ W2 = G" + n(F - G)G'''^ + ^{n\F - Gf - r2)G''~\ 
(b) In the important case G = F, the formulas above become somewhat simpler. Here, we derive 

V, = 0, V2 = -^T2, Vs = -^Ts, y4 = ir2-ir4, (ii) 
V5 = -r2r3--r5, = -— + -r2r4 + — - -Fs, (12) 

5 g i d g &, b 48 2 ^ 8 18 ^ 6 ' ^ ^ 

V7 = -^riF3 + ^r2r5 + ^F3F4-^F7, (13) 

Vs = —ri - — F^F4 - — FsF^ + — F2F6 + — F3F5 + — F? - -Fg, (14) 
^ 384 2 32 2 36 ^ 12 " 15 32 ^ 8 ^ ^ 

which, in particular, leads to Wq = Wi = F", 

W2=r'- \t2F^~\ Ws = r'- \t2F^~^ - ^r3F"~'. 

Letting £ = 1 and a = 0, we obtain under the present assumption that 

n^^--^" ' if ,i<(2eci)-\ (15) 



where 



m = max —r— (16) 

fcen\i k 



(see comment after Theorem 12. ip can be estimated with (jl8p below. 
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(c) Let us assume that, for each j G n, Fj <C G and let fj be a G-density of Fj. Set / = ^ Sj=i 
If in Theorem I2.ir b) we choose suitable Bq, Bi, . . . , Bn G {0,3C}, it then follows that, for 
i G Uq, 



L„„„{2„£(7-lfdG,^||F-Gf} 
+ E->""{^ £(/j-l)'<iG,^||F,-G|r}. (17) 

Prom the definition of r/^ it is clear that, if G = Fi = ■ ■ ■ = Fn, then rj^ = for each i € Uq- 
The inequalities (jl7p and (jlOp reflect this fact. Moreover, in view of these bounds, if G ~ F in 
some sense and if the Fi , . . . , Fn, are not to o different, then a large n leads to a small bound. 



Speaking in terms of lBarbour et al\ (|l992l . Introduction), our bound contains a magic factor 
(cf. Section [L3] above) . 

(d) If G = F, then, for each j G n, we clearly have Fj <C G and therefore a G-density fj of Fj 
exists. In this case, pT|) reduces to 



^ j2'^in[^ju3 - 1)' dF, ^||F,- - F||'}, {I G n). (18) 
We note that, in (fTSj) . /^(/j — 1)^ dF is finite for all j G n, which follows from 



/ /|dF= / /,dF,-^n / /,dF = n. 
JX Jx JX 



One might ask whether the singularity in the right-hand side of Q can be removed. The 
following theorem shows, that this is possible, if we enlarge the leading absolute constant and 
replace r/^ with rj^ (or with 7]^ in the case G = F). 

Theorem 2.2 Let the notation of Theorem \2.1\ he valid. 

(a) Let i £ Hq and let ug G (0, oo) be the smallest possible constant such that, without any 
restriction on r/g, 

n 

\\llFj-Wi\\^ue,r^/\ (19) 

i=i 

We have 

(2eci)(^+i)/2 ^ ^ 

Ui ^ ^-^ , 20 

l-xi 

where X£ G (0, 1) is the unique positive solution of the equation x^'^^ + x/2 = 1. By (|20p . we 
get uq ^ 5.9, ui ^ 17.3, U2 ^ 44.5, and u^, ^ 107.5. 

(b) Let £ £ n and let U£ G (0, oo) be the smallest possible constant such that, under the assumption 
G = F and without any restriction on rj^ , 

n 

\\j\F,-w\\^u,r,^^^'^/\ (21) 
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Then we get 



5i ^ 10.94, ila ^ 31.5, ^ 82.2, ue , Uen\3), 

I- Xi 

where X£ G (0, 1) is the unique positive solution of the equation x^'^^ — x^/2 + x = 1. 



(22) 



Remark 2.2 (a) If 779, resp. ry^^, is sufficiently small, the bounds given in Theorem 12.21 can be 
further improved as follows from Theorem 12.11 and Lemma 14.51 below. In particular, in the 
case G = F, we have (cf. proof of Theorem 12. 2p 



(23) 



n 

i=i 


-y 




^-2,2 + 


n 

|n^> 


-W2 


^ (1 + n2^)r/i. 


n 

II n^'. 


-W2 




^-3,3 + 


n 

|n^. 

i=i 




^ (\/3 + U3^)rjl^'^ 


n 

II n^. 


-W3 




1 2 1 

8-4,2 + 1 


n 

|n^. 


-W4 


^ (2 + U4^)r/^ 



In view of (j23|) . one may conjecture that ui ^ 1. Indeed, this is correct and follows from the 
simple observation that, for j£ = Z, n G 2N, Fi = ■ ■ ■ = F„/2 = -^O) Fn/2+1 = ' " = Fn = h, 
we have 

n \ 1 
,n/2y 2", 



2 1 



2, (n ^ cxd) 



and, by ^2) and (US}, || n"=i -P'i - ^ ui 7?i ^ 2ui. 
(b) From dUD and it follows that 



in^.-^' 



^ 2ni max 



lYdF. 



(c) It is unclear, whether it is possible to remove the singularity in ([SD for any a > 0. Indeed, 
since the denominator of the right-hand side of ([8]) contains rji and not rji^a, we cannot argue 
as in the proof of Theorem [ 



Example 2.1 In the situation of Theorem l2.1l let us assume that Fj = Ylt=oPj,rHr, (j S n, d G N) 
and G = F = Ylt=oPr^r, where Hq, . . . , Hd e and for r G d^, pj^r £ [0, 1] with Ylt=oPj:r ~ ^ 
and Pr = ^ Yl']=iPj,r > 0. Then, for each r & d^, Hr has a F-density hr and we may assume that 
Ylt=oPrhr = 1- Consequently, Fj has the F-density fj := Ylt=oPj,r^r, (j G r)- Using the simple 
inequality 

\2 '^2 

" h). (24) 



^fe^^E^' (a. G [0,00), (0,00) for re do 



r=0 
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we obtain, for j S n, 

ifj-lfdF = f ffdF-l= f (^r=0Pj,rhr)\ -p_^ 
X Jx Jx Y.r=QPr^r 

" ^ hK>^} ^ Pr- Pr ' 

Further, we have 

d 

W^j ~ \Pr-Pj,r\- 
r=0 

Therefore, in this context, (|18p imphes that, for i £ n, 



Using (pH) . we get 



j=l r=0 r=0 



{Po-Pjflf _ i'^ - Po){J2t=liPr -Pj,r)? ^ /' 1 V" iPr-Pj,rf 



PO PoY.t=lPr ^ ' ^1 Pr 

and hence 

We note that (j26p is non-trivial in the sense that the sum on the right-hand side does not contain 
the summand for r = 0. In view of (j2ip . (j22p . and (j25p with ^ = 1, and (j26p . we see that, in ([7]), 
the factor Csd can be replaced with 2ui, which in turn is bounded by 21.88. We note that, if the 
Hr are given as in ([1]) then (j25p and (|18p coincide. But if //q ~ • • • ~ in some sense, then (j25p 
can be much worse than p8p and should therefore not be used in general. 

The next proposition shows that, as claimed above, sometimes r]^^, (a > 0) has a better order 
than r]^. Here, we consider the case of symmetric distributions Fi, . . . , F„ G ^ with finite support. 
For simplicity, we assume that G = F. 

Proposition 2.1 Let the notation from Theorem \2.1\ hold. Further, let b G N, xi, . . . , € X\ {0}, 
= Pjfil + Ylr=iPjA^-xr + Ixr) e ^, U G nd, and G = F = PqI + J2l=iPr{I~xr + Ixr), where 
Pj^r G [0, 1] with pjfi + 2 Ylr=i Pj,r = 1 '^'^^ Pr = ^ Sj=i ^'j'^ ^ ^' ^ ^o) ■ ^'^'^ ^ ^ R; we then have 

V,, ^ t min f'^^' + 2 V - + i: " f ) , 

^PO ^ M ^ Pr ^ 

(£ + 1)2 (l^O ~ Pj.ol +'^Y.^Pr- PJ^rl) }• (27) 



OE 



r=l 

,-2 



We note that, in contrast to (I18p . the bound in (I27p has the better magic factor n . Hence, in the 
situation of Proposition 12.11 estimate ([8]) with a = 1 should be preferred over dH). 
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3 Numerical examples 

In what follows, we compare the available bounds in the multinomial approximation of the gen- 
eralized multinomial distribution. We assume the notation given in ([T|) with d = 10. Further let 
i = 1. The following two examples show that the results of the present paper can be considerably 
sharper than the bounds from the literature discussed in Section 11.31 

Example 3.1 For j S n, let pj^r = (r)^^(l ~ ^ do) be the binomial counting density 

with number of trials d and success probability qj = 0.4 + jj^^, where a ^ 1. Clearly we 
have Qj G (0.4, 0.5] for all j £ n. We emphasize that, with this definition, Fj is not a binomial 
distribution. Further, if a or n is large, then pj^r should be close to Pj. for a sufficient number of 
j £ d and r £ d^, so that we expect a small distance || 11^=1 ~ -^11 here. This is reflected in the 
bounds, given in Table [TJ 



Table 1: Numerical bounds for the distance in Example 13.11 



n 


a 


Ci 


C2 


dSD 


© 


dSD 


(|23]) fc ([25]) 


([T5]) fc ([25]) 


100 


1 


111.4 


15590.9 


n.a. 


> 2 


n.a. 


0.197438 


0.173503 


1000 


1 


145.7 


26444.8 


n.a. 


> 2 


n.a. 


0.026902 


0.032981 


100 


2 


154.6 


29809.2 


n.a. 


0.107737 


0.034777 


0.000366 


0.000954 


1000 


2 


156.3 


30455.0 


n.a. 


0.110925 


0.035914 


0.000037 


0.000120 



Note that the bounds for the distance are always rounded up. Further, as the distance is always 
bounded by 2, larger bounds are omitted. The entry "n.a." means "not available" and describes 
a situation, where the bound cannot be used, since the respective condition does not hold. In all 
cases, the quantities Ci and C2 (see (|4|) for the definition) are quite large, which explains that the 
condition for ([3|) is not valid here. This is due to the fact that, in each case, some of the (r G dg) 
are quite small. E.g. see Table [2] for the case n = 100 and a = 1. 





Table 2: 


Point probabilities of F in 


Example 13.11 when n 


= 100, a = 1 




r 





1 


2 


3 


4 


5 


Pr 


0.00416 


0.03012 


0.09851 


0.19175 


0.24611 


0.21781 


r 


6 


7 


8 


9 


10 




Pr 


0.13473 


0.05757 


0.01628 


0.00276 


0.00021 





In the next example, we discuss a situation, where ([3|) gives non-trivial bounds. 

Example 3.2 For j £ n and r £ d^, let 

_ l + + r)/{b{n + d)) 

^'''~Ei=o(i + 0- + n)/(Kn + rf)))' 

where 6^1. Similarly as in Example 13.11 for large n or b, we expect good approximation, which 
indeed is refiected in the bounds for || YVj=i ~ ^|| given in Table [3l 
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Table 3: Numerical bounds for the distance in Example 13.21 



n 


b 


dSD 


dSD 


dSD 


dlSD & ([25]) 


(dH) & ([25]) 


100 


1 


0.325253 


0.008310 


0.002337 


0.000030 


0.000098 


1000 


1 


0.118021 


0.000119 


0.000033 


3.9 X 10^^ 


1.5 X 10-^ 


100 


2 


0.112763 


0.000978 


0.000267 


3.3 X 10~^ 


1.2 X 10"^ 


1000 


2 


0.040581 


0.000014 


3.8 X 10"*^ 


4.4 X 10"^ 


1.7 X lO"'^ 



In contrast to Example 13.11 in each case the values p^, (r G dg) are quite similar, which implies 
that the condition for Q is valid. E.g. see Table H] for the case n = 100 and 6 = 1. 



Table 4: Point probabilities of F in Example 13.21 when n = 100, 6=1 



r 





1 


2 


3 


4 


5 


Pr 


0.08807 


0.08864 


0.08921 


0.08978 


0.09034 


0.09091 


r 


6 


7 


8 


9 


10 




Pr 


0.09148 


0.09204 


0.09261 


0.09318 


0.09374 





In what follows, we discuss an example, where the distance can actually be evaluated. 

Example 3.3 Suppose now that, in Example 13. H we change the measures Hr to = on M for 
r G df), i.e. all distributions Fi, . . . , Fn, F are one-dimensional. Then, using a computer, it is not 
difficult to get the exact numerical value for the distance, see Table [5l 



Table 5: Exact numerical values for the distance (cf. with Table [T]) 



n 


100 


1000 


100 


1000 


a 


1 


1 


2 


2 


lin-=ii^.-^ll 


0.007152 


0.001653 


5.9 X 10"^ 


7.6 X 10"^ 



A basic property of the total variation distance tells us that, for distributions Hr £ J^, (r £ ([q) in 
the case of a general measurable Abelian group, we have 

d n d d 



n 



d 

E 

r=0 



Pj,rHr 



r=0 



I + ^Prh 



r=l 



r=l 



(28) 



This can easily be seen by writing the difference of the measures on the left-hand side as a polynomial 
in Hr, G dg) then applying the triangle inequality. As a consequence of (128]) . each bound 
from Table [1] is valid here as well. A comparison shows that the bounds are getting closer to the 
actual distance as n or a is becoming large. For example, the bounds from (j23p & ()25p are about 
27.6, 16.3, 6.2, and 4.9 times higher, respectively, than the values from Table [5j 

We can apply this idea to Example 13.21 as well: if we again change the measures H^ to Hj. = 1^ 
for r G get the exact values of Table El A comparison with Table [3] shows that the bounds 

from ([23]) & (f25ll are about 4.8 to 4.0 times higher than these values. 
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Table 6: Exact numerical values for the distance (cf. with Table [3]) 



n 


100 


1000 


100 


1000 


b 


1 


1 


2 


2 


lin-=ii^.-^ll 


6.3 X lO"*' 


9.1 X 10-8 


7.4 X 10"^ 


1.1 X 10-s 



4 Proofs 

4.1 Smoothness estimates for convolutions 

In what follows, we use the standard multi-index notation: For z = {zi, . . . , Zd) G C^, {d G N) and 



w 



{wi, . . . , Wd) € Zl, we set 



Ylr=i^r, and ty! = nr=i''^f'- Similarly, for 



V = {Vi,..., Vd) G M'^, set y"' = nr=i K"""- For u,?/; e Z^, we write u ^ w in the case that 
Vr ^ tL'r for all r £ d; let V A w = {vi A wi, . . . , Vd A Wd)- Sums over v, v, and w are takeii over 
subsets of zi as indicated. The following lemma is a counterpart of Lemma 5 in iRood ( 200ll ). 



Lemma 4.1 Let k,n Z+, d € N, and £ M. for v € Z'J, wii/i |f| = A;. Let X = {Xr)red be 
random vector in with E[(^^^-^ < oo and put Xq = Ylt=i -^r- Let p = {pr)red G (0, 1) 

such that po = l- Yfr=iPr e (0) !)• Further, let H = {Hr)r&d G Hq G J^, 



a a a j 

G = Y,PrHr eJ^, Ul=Y,^ UiHr - HoY^ , ^2 = e( Xr{Hr - Ho)) 



T 

r=0 |i;|=fc r=l r=l 

where, in the definition of U2, the expectation is defined setwise. Then we have 

'I / w\{k-\w\)\[ ^ a„ A ^Vr\^^\^/'^ 



|t)|=fc r=l 



n + k 
k 



r=0 



(29) 

(30) 



where the random vector Y 
Proof. Let 

Mult(u), n,p) 



{Yr)red is an independent copy of X and Yq = ^^^^ 



\w\ ^ n. 



w\in^\w\y. P Po ' IX uy c 

0, otherwise 

denote the multinomial counting density with parameters n and p. For / : Z'^ — > M and r £ d, let 
Ar/ : Z*^ — > M with (Ar/)(ti;) = f {w — e,^) — f {w) for w S Z*^. Products and powers of A-operators 
are understood in the sense of composition. Further, let A^/ = /. Clearly, Aj.^ Ar2/ = Aj-jA^^/ for 
ri,r2 G d. For t; G Z^, let A''/ = A^i • • • A^V- We set A''Mult(i«, n,p) = ( A^Mul t (-, n, for 
u; G Z^. We use the following properties of the multinomial distribution (see IRoosI (|2001| . formulas 
(20), (21), and (4))): For v G Z^, 

d 

A''Mult{w,n,p)H'" 

|ui|^n+|?j| r=l 



A^Multiw, n,p) H"" H'^+\''\-\'"\ = JJ(F, - HqY^ (31) 
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and, for v,w £ 



A''Mult(w, n,p) = Kraw(f ; n + |z;|,p)Mult(i(;, n + \v\,p) 



v\ n! 



{n + \v\)\p''p^^^ 



(32) 



where 



Kr:aw(v;w,n,p) = 



n — \w\ \ \v — v\\ {—p)^ ^Pq''' 
Iv — f I / — w)! 



n 

r=l 



(33) 



is a Krawtchouk polynomial of degree v. Note that there is another set of Krawtchouk polynomials 
which forms, together with the one fror n (1331). a bi-ort hogonal system of polynomials with re spect 



to the multinomial distribution (see also iTratnikI (jl989l )). Prom the more general Lemma 2 in 
(|200lh . it follows that, ioi v,v with | v\ = \v\, we have 



Roos 



Mult(w, n + \v\,p) Kraw(v; w,n + \v\,p) Kraw(f ; w,n + \v\,p) 

\w\!^n+\v\ 

(n + |v|)! |v - u'|!p^+'^-"'p[)'"+''' 



wl n\ (v — w)l {v — w)\ 



(34) 



We note that the right-hand side of (j34p is always positive, which shows that, if d ^ 2, then the 
Krawtchouk polynomials given above are not orthogonal with respect to the multinomial distri- 
bution. However, we do not need such a property. Using ([3T]) . ([32]) . Cauchy's inequality, we now 
obtain 



= \\Y.^ Yl A^Mult(u;,n,p)F^F, 



V. 

\v\=k |ui|^7i+|d| 







n! 



{n + k)] 



{n + k)] 



Mult{w, n + k,p)\ ^ ^ ^ Kraw(t;; w,n + k, p) 



\v\=k 



pup^ 



Mult(w, n + k,p)\ — ^Kraw(?;; n -|- A;,p) 



2\ 1/2 



T. 



ui6 



Using (jHlj) . we get 



n! 



(n A;) 



\v\=k\v\=k ui^dAi? 



E 



(n + A:)! (A: - |tf|)!p^+''-"'p|7'+''\i/2 



E w!(fc - |w|)! r ^ TT / ^'r \ TT 

|i,|=fc r=l ^ '^^ |5J|=fc 

w!(A; — |w|)! r ■^-^ a,, -A- / \ 1 2\ 1/2 



w\ n\ {v — w)\{v — w)\ 
d 



Jin + k)\ ^ n«'nn~''"' f! -^-^ ytfr/ f! -^-^ \ Wr 

V ^ ^ |ii)|^fc -f^ -f^O \v\=k r=l ^ ' \v\=k r=\ 



^(^^Tfc)! 



1/2 



E wii^ft; — \u}\)'. r a„ T-r 
|«)|^fc /^O |i)|=fc r-=l 



Vr 
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Inequality (l29|) is shown. Since U2 = E|d=fc I Ellil^r " ^o)^^ (EH) gives 



n\ k\ 



E 



w\{k — \w\)\ 



k-\w\ y\ W 

\w\^k P PO \v\=k r=l 



2\ 1/2 



For ^ A;, we have 

\v\=k r=l 



u)|^fc P |»;|=fc r=l ^ \v\=k r=l ^ 



1/2 



X'"X, 



A;— 




wl {k — \w\)l 



Indeed, this follows from the identity theorem for power series taking into account the fohowing 
equality of the corresponding generating functions 



k=0 \v\=k ' r=l 

From the above, we get 

/ I ; \ -1/5 
n + k\ ' 



■Wl 



r VW V 



k—\w\ 



fc=|ui| 



w\ {k — \w\)l 



{z £ C). 



k 



n + k 
k 



E 



-1/2 



d 



k\ 



XQYQ\k-\M 



EE 



w\ {k — \w\)l \ po 



n 

r=l 



Pr 



1/2 



r=0 



Pr 



which completes the proof of ([30 



□ 



T he foh owing lemma is an important application of Lemma |4. II and generahzes formula (37) in 



RoosI (|200d ). Another application is given in the proof of Proposition I2.H see Section [4.31 below. 



Lemma 4.2 Let € N, n G G G JT, and U G Ai, where we assume that \U\ <^ G and 
that U{X) = 0; let denote any Radon- Nikodym densities of with respect to G and put 
f = f+-f-. Then 

k/2 



n + k^ "^^^ 
k 



rdG 



(35) 



Proof. If j f^dG = 00, then (j35p is trivial. In what follows, we assume that / /2 dG < 00. Let 
e G (0, 1) be fixed. Then af^^ G [0, cxd), {j G Z+) and pairwise disjoint Bj^ir G A, (j G Z+) exist such 
that 

00 00 

|Ji3,-, = X, :=J24e^iBj,e), and ^ /± - ^ e. 

j=0 j=0 

Here l(^) is the indicator function of a set A. Let be the measures on {X,A) with G-densities 
ff. This implies that Uf = Yl'jLo Oife Hj,e, where, for j G Z^., 



Qte = 4sG{B,,e) and i/,,. 



G(B,,,n ■)/G{Bj,e), if G(i?,, ,) > 0, 
I, otherwise. 
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Set = qj~^ ~ ~ ft ~ fs ' ~ ~ ■ r^oie that the latter equaUty 

indeed indicates the Hahn- Jordan decomposition of Ue. Then ||f/e|| ^ and ||?7 — C/gH = 

/ ((/+ - ft) + if- - fr)) dG ^ e, giving 



Hence 



Since //dG = t/(X)=0, we have ET=o'l,,e = /(/+-/-)dG 
for each m G N, 

k / ™ ^k 



^ lie/ — C/eII ^ e, and therefore, 



I ( Ij'^^i'' 



j=0 

Hence, we obtain 

oo 



j=0 j=m+l j=0 



k-l 



j=0 j=0 

oo oo 

^k[e + 2 |g,,,|)(2^|g,-,| 

j=m+l j=0 



k-l 



k 



j=l j=0 

From ()30p . it follows that the norm term on the right-hand side is bounded from above by 



n + k\ "^^^ 
k 



jemo:G(Bj,e)>0 



G{B 



3,e) 



n + k 
k 

n + k 



j=0 



Letting m ^ oo, we obtain 



oo J 1 / r 



3=0 



-1/2 



/'dG 



k/2 



Since 



I if - ft) dG| ^ / 1/ - M(l/I + I/el) dG ^ e{\\U\\ + \\Ue\\) ^ 2e\\U\ 
we obtain (135]) by letting e — > 0. This completes the proof. 



□ 



It may happen that the assumption in Lemma 14.21 does not hold directly. However, this can 
sometimes be overcome by shifting U. The following corollary is needed in the proof of Theorem 1 2. 11 



Corollary 4.1 Let n G Z+, G £ J^, Ui,U2 e M, and U = U1 + U2. We assume that \U2\ <C G 
and that both / 0. Put = U2/\\U2\\- Let denote any Radon- Nikodym densities of 
with respect to G and set f = f^ — f" ■ Then 

\\UG-\\ ^ ||f/i|| + |C/2(X)| + ^^£i^M( / fdG)'". (36) 
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Proof. The assertion easily follows from the triangle inequality, Lemma 14.21 ^-iid the simple fact 
that U2 = [/2(^)f^2 + (ll^^ll ^ ||^2~||)(^2' ~ f^2~)' where r denotes + or — according to whether 

II ^2" II > 11^2^11 °^ '^O*- ^ 



Remark 4.1 (a) Let the assumptions of Corollary 14.11 hold. If is a cj-finite measure on X and 
if G <C /i, then G and have ;U-densities v and g^, say, and, letting g = — g~ , we can 
write lfdG = f^^^^^ g'^y-^ dfi. 

(b) Sometimes it is useful to simplify further the bound ()36p by using the following inequality: 
(11^2^11 ||f^2 ||)^ / ^ jh^dG, where h = - h~ and denote any G-densities 
of U2 ■ Indeed, this follows from the representation 

J Ja IIC/2 II •^^'^ 11^2 II 

whenever A with U2{A) = U2[A'') = 0. 

The next corollary is an extension of Lemma 14.21 to compound distributions and may be par- 
ticularly useful in the compound Poisson approximation. 



Corollary 4.2 Let /c G N, G G JT, and U G where we assume that \ U\ <^ G and that U{X) = 0; 
let denote any Radon- Nikodym densities of with respect to G and put f = f~^ — f~- Let 
N be a random variable in Z+ and (p{z) = E[z^], {z G C, \z\ ^ 1) be its generating function. Set 
ip{G) = E[G^] G J^, where the expectation is defined setwise. Then we have 



1/2 



f^dG 



k/2 



\\u''(f{G)\\ ^ {k J x'^- V(l - a;) dx) 
If N has Poisson distribution exp(t(/i — /)) with t G (0, 00), then 

P'^iG)\\ ^-L,/klp(NTk){ I f'dG) 
Proof. Using the triangle inequality. Lemma 14.21 and Jensen's inequality, we obtain 



(37) 



k/2 



(38) 



|[/V(G)|| ^ E||?7'=G^|| ^ (E 



N + k 
k 



-1, 



1/2 



f^dG 



k/2 



The integral representation of the beta function implies that E(^^'^) = k x'^^^ip^l — x)dx, 
which, in turn, leads to (j37p . Inequality (j38p easily follows from (j37p and the series representation 
of the lower incomplete gamma function. □ 
We note that (I38|) is comparable to previous results of iRoosI (|2003l . Lemma 2) but is however 



much better because of the more general assumptions used in Corollary 14.21 
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4.2 A general lemma 

The results of Section [2] are based on the following general lemma. Here, a distribution G £ is 
called infinitely divisible if, for each n G N, there exists a Gn £ such that GJ^ = G. We note 
that, in general, such a n-th root G„ need not be unique (see iHeveii (|l977l . proof of Theorem 3.5.15, 
pp. 222-223)); let G^/" denote any fixed n-th root of G. 

Lemma 4.3 Let n e N, Fi, . . . ,Fn,G e J", Li, . . . ,Ln e M. SetL = ^ ^"^^ Lj, Kj = FjQ-^i , 
(j ^n),Ko = Ge-^, 



JCn: \J\=kjeJ k=0 

f (i^j - Ko)(G"-'=)l/'=e^ if G is infinitely divisible, , , 



otherwise, 



^Mj-fc , (A; G n,m G N). 
Lei he defined as in Theorem \2.1[ Then, for i G Uq, 



\llF,-We 



A 



The following two examples show possible applications of Lemma 14.31 As a byproduct, results 
in the compound Poisson approximations can be derived. 

Example 4.1 If we consider the case Li = • • • = L„ = 0, we see that Theorem 12. If a) is a direct 
consequence of Lemma 14.31 (cf . proof of Theorem I2.ip . 

Example 4.2 Suppose that, for j G n, Hj G J^, pj G [0, 1], Lj = Pj{Hj — /), and Fj = I + Lj. Put 
L = Sj=i ^^'-^ G = e^. Then Lemma 14.31 implies that 



n n n , 

|n^.-G"||«E^(EiK4)- 

j=i k=i ■ j=i 



(39) 



where, for j, k £ n, 



{I + p^{H,-I))e 
iKj-Ko){G 



-Pj{H,-I) 



Kn = L 



"n—k\l/kL 



((/ + L,)e- 



I) exp ( -L 



n- 



In principle, (j39p is the same as estimate (26) in iRood (120031) . The approach used there is based 
on a slight modification of an expansion due to iKerstanI (jl964l ) . It is however not sufficient to get 
the results of the present paper. 

For the proof of Lemma 14.31 we use formal power series over M . In the following lemma, some 
basic properties in connection with the norm on M are summarized. The proof is simple and 
therefore omitted. 
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Lemma 4.4 For n € N and k £ n, let '4)^^\z) = Yl!f=Q^j,kZ^ , {Wj.k S M) he a formal power 
series over J\A with variable z and let CoeS{z^ ,'ip^^\z)) beitsjth coefficient Wj^k- Further, consider 
the formal power series ip^j^\z) = l^^o ll^i.^ll •^'^ indijj^^\z) = Yl'jLoO'j,k z^ foraj^k £ oo) 
and k € n. Then, for j £ , 



\CoeS{z^,ijf\z))\\ = CoeS{z^,i;\^\z)), 



n 



Coeff ( n (^) ) ^ Coeff (z^ ,ll^Pl}\z)) ^ Coeff ( z^ J] vf ^ (z) 



^ Coeff (z^ J] 

fc=l k=l k=l 

Proof of Lemma 14.31 We first note that 



i=i i=l fc=0 A:=0 

For k n, let A(n, /c) = or A(n, A;) = n — k — k[{n — k)/k\ according to whether G is infinitely 
divisible or not. For i €z Hq, we obtain 

ran n 

j=l k=l+l k=l+lJQn:.\J\=k j£J 

n 

= Coeff(z^^fc(z))G^("■'=), (40) 

k=e+i 

where ipk{z) = 0^=1 (-^ + ^j,kz) is regarded as a formal power series for k £ n. It should be 
mentioned that it is essential here to extract the kth. coefficient of a formal power series which itself 
depends on k. By Lemma 14.41 ^oi k £ n, we get 

n j.k 



|Coeff(/,?/^fc(z))|| ^ Coeff (^z\JJ(l + ||Mj-fc||z)j ^ Coeff (z^ e'^^'i^) = (41) 



On the other hand, using 



Mz) = exp(YMj,kz)ll{e-^^''''iI + M^^kz)) 
i=i i=i 



n oo , 

1 — m 



exp (^M,,,.) n [ E 



j=l j=l m=0 

we derive 

n 

||Coeff(z^Vfc(^))|K Coeff (z^e^'=^^^^(|l^^>ll^))' (^2) 

where, for y £ C, 



oo 



= y ^^y- = 2-(l-y)e^ = l + ^ + ^ + ^ + .... 
^^^^ ^ m! ^ V 2 3 8 

From the definition of Ci, we obtain that 

\giy)\<gi\y\)^e''\y\\ (43) 
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Here, we note that h{x) := ln(2— (1— x)e^')/x^ for x £ (0, oo) attains its maximum = 0.694025 
at point xq = 0.936219 .... This can easily be shown using the representation 



h{x) = \ ^ ln(2 - (1 - tx)e*^) dt = C — — ^- dt, 

^ ' x^ Jq dt ^ ^ ' ' Jq 2e-*^ -1 + tx ' 

which, after differentiation of the integrand, leads to a useful integral formula of the derivative 

, 1 r i^(2e-*-l) , 

h'{x) = ^ —\ ^dt. 

^ ^ x^ Jo (2e-* - 1 + t)2 

As a consequence, we learn that h'{x) = has exactly one positive solution x = xq, which can 
easily be calculated numerically. Let 

Bessel(0; y) = ^ Tl)^ " 2^ / """"P^^ """"^^^^ ^ 

be the modified Bessel function of first kind and order 0. Using (j42p . Cauchy's integral formula, 
and (03]), we derive, for /c € n and arbitrary E (0, oo), 

||Coeff(z^V^fc(z))|| ^ r e-i'=*(j]5(||M,,fe||iifce'*))exp(i7fciifce'*)dt 



— ^ / exp(i7fciifcCos(t))dt ]^ c/(||Af,-fc||i?fe) 



^ 7^Bessel(0;i:^fci?A:)exp(cii/fe2^fc) 



;2 



where (p{x) = Bessel(0;x)e ^'^/^ ^ 1, (x G M). Choosing 
we get 



|Coeff(z^V^,(z))|| ^ (^(^ + ^,2))'^'- (44) 



Taking into account (jlDj) . the fact that A(n, A;) G ngj well as ([^T]) and (jM]), we obtain 
InF^-H-.l < Y. ||Coeff(.'.V,.(.))||< E [(^(^7 + -'.^)) '^^ 

j=l k=t+l k=i+l 1 

The proof is completed. □ 
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4.3 Remaining proofs 

Proof of Theorem 12.11 Part (a) follows from Lemma |4.3[ Indeed, for tj^ < (2ec^)~^, we have 



E 

k=i+l 



k \4c 



k=e+i 



(l-v/2^)/3+i- 

Here we used that, for x G [0, 1), 



Part (b) is shown by using Corollary 14. II together with Remark 14. 1[ In fact, for j, E n, we obtain 



l/2\ 2 



since, for k £ n, 



Similarly, we have 



n — k 
k 



^ max{n — k,k} ^ n 



i:;| = ||X^M,-fc||'^ (n||(F-G)|Bj|+n|(F-G)(Bo)| + V2fc^( / (/o-l)'dG 



l/2\ 2 



This yields (fTOj) and completes the proof. 

For the proof of Theorem 12. 2| we need the following lemma. 



□ 



Lemma 4.5 Let n £ N, Li, . . . , Ln £ Ai with Yll=i — ^' '^''^d, for k £ and m € N, 



Then we have 



JCn: \J\=kj£J j=l 



\\V2\\^l^2, \\V3\\^l^3, IM^l^, 

zoo 

\\y4^l^2^3, iM^r^^l ll^riK^^i^s. 



144 



24 
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Proof. The first two inequalities are easy. Taking into account pT]) - p^ . it is not difficult to show 
that, letting = Ej=i(--^i)"', ("^ ^ N), 



- [r^-r4]-r4 



^ o(l|ri-r4|| + ||r4||) 



[rsFs - r5] - -Ts 

5 







1 



||3[f 3 - 3r2r4 + 2r6] - 9[r2f 4 - Fg] - 8[r2 - Fe] + rg] 

^ {Ml + 8[7?i - M - 2i?6) ^ 



144 
1 

840 
1 



144 



'2' 



35[r^r3 - 2r2r5 - r3r4 + 2r7] - I4[r2r5 - Tr] - 35[r3r4 - Tr] + Tt 



840 
1 



= ^ l^8[^?ii93 - 2?92i95 + ^97] + - 6^97) ^ ^i92^?3. 

Observe that, in order to obtain good constants, a convenient grouping of terms is essential. Further, 
for the bound of 11 Veil, we used the inequality {'d'l^ — 'dQ)/'d'2 ^ 4~^, which can be proved by using 



^3 



n 2 " " 

(X]^''^) ^ ^x]{l-Xj) =: gn{ixi,...,Xn)), 



where Xj = ||Lj|| ) , and the fact that the functions forn E {3, 4, ... } 

and {xi, . . . , Xn) G [0, 1]" with Yl^=i — ^ satisfy 

9n{{xi,...,Xn)) ^ gn-l{{xi + X2, X3, . . . , rE„)), 

whenever ^ xi ^ . . . ^ ^ 1. This completes the proof of the lemma. □ 

Proof of Theorem 12.21 In order to prove the assertions, we need a further bound. In fact, 
similarly as in the proof of Lemma 14.31 we get that, for i £ Uq, 



\l{Fj-We < 1 + IIVF^II ^ 2+ ||^Coeff(^^^]J (/ + Mj-fcz)jG^("''^) 
j=i k=i j=i 



2 + 



k=l 



k=l 

2-t- t'+^ 
l-t ' 



where t = y^2e ?7g. Similarly, \i G = F, then, for i £ n, 

2-2?+?-?+! 



\llF,-We ^2 + J2t 



7k 



k=2 



l-t 



where t = ^J2e rj^. We now prove (a). Let i £ tiq- If t £ [0,X£], then ([9]) gives 



\llF,-W, 



1 — t 1 — X£ 



22 



On the other hand, if t £ {x£, oo), then 



\llF,-W, 



1 - t 



1 — Xf\ 



xe\t'+\l-t) J\x'+\l-xe) ^ l-xe' 



since = t^^^^ + Z^jii * is decreasing on (0, oo). This yields (fT9]) and (f20]) . The proof of 

(a) is easily completed. Let us now show (b). Set G = F. Similarly to the above, one can show 
that, for i G n, 



\YlF,-We 



I- Xi 



This proves one part of ()22p . Using the norm inequalities in Lemma 14.51 and Q, we derive, for 

ie3, 



i=i i=i 
where, for x S [0, {2eci)~^), 



+ 



Wj-We 



j=i k=e+i 



Ci{x) = x + VSx^^^ + 2x^ + ^x^/^ + ^x^ + ^ 

6 2 24 



7/2 (2eci x) 



l-y2i 



6 2 24 l-y/2i^' 



w N n 2 55/2 ,2 15 3 77/2 ,2 2ecix^ 

Cafx) = 2x^H x^'^ -\ x^H x^'^ + — } ' 

' 6 2 24 l-y2i^ 

Note that, for ^ € 3, we have ^(^i) ^ 2-2?+?-?+^ ^ •£ g^j^^-^ ^^^jy -f ^ [0,5^], where si = 
0.182839 . . . , S2 = 0.196439 . . . , and S3 = 0.205094 .... If ?7i G [0, s^], then || n"=i - ^ 
CK^i) ^ If '?i £ (s^^oo), then, letting ti = y^2eqs£, 



i=i 



^ ^+1 ^ 2t + P t^+-^ ^ 2 2ti + tj 



(<?+l)/2 QjSi) 



.(^+l)/2 ■ 



Numerical calculations give the bounds for U£, {£ £ 3) as claimed in ([22l) . This completes the 
proof. □ 



Proof of Proposition 12.11 Consider fixed j, ken. Let d = 2b, p = {pi, ... ,Ph,pi, ...,%) G Z'^, 
and p = [(n — k)/k\. Further, for v G Z'|_ with |f| = 2, let a^, = p^ ~ Pj,r if = I'b+r = 1 and 
ttj] = otherwise. Let 

j r £b, 



Hn=L 



Hr 



•^r — b ' 



r e d\b. 



Closeness of convolutions 
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Then we have F = Ylt=oPr^r and 



iPj,0 - Po)I + ^{Pj,r - Pr){l-xr + Ixr) 
r=l 

b b 
^{Pj,r - PrW^Xr + Ixr - 2/) = ^{Pr " Pj,r){I-Xr " " I) 



r=l 



\v\=2 ' r=l 

Here and henceforth, sums over v and w are taken over subsets of as indicated. In particular, 
we obtain 

6 



= WiFj - F)F'\\^ ^ \\F, - Ff ^ {\Pj,o-Po\ + 2^ 
On the other hand, in view of 



Pr 



r=l 



\Mj,k\f = \\iF,-F)F'\f 



E^^U(Hr-H.r){EPrHr 



\v\=2 r=l 



r=0 



we see that ([22D can be apphed, which together with the simple fact that 



(p+2)! ^ (p+1) 



gives 



M 



i,*: I 



E 



w\{2 - \w\)l r 



(p + 2)! ^ »5i«:f52-|«,| 
^'^ ' \w\<^2 P PO \v\=2 r=l 



[Eirn 



< / 2 r g^,! 2 1 [ ay -At f Vr\ 

1^1=2 \w\=l^ \v\=2 r=l ^ 



+ 



|to|=2-^ \v\=2 ' r=l 



The special definition of a^,, (v G Z^, |v| = 2) implies that = for r(l),r(2) G d with 

|r(l) — r(2)| 7^ 6 and therefore the terms on the right-hand side can be evaluated as follows: 

b „ , _ .2 



\v\=2 



r=l 



(Pj.o -P oY 

'^pI 



|«;|=1 



\v\=2 r=l 



Ew^ r v;^ TT- 
■pui [ Z-/ til J- J- 



Wr 

Vr 
Wr 



V— — [ V )] = 2 V 



1^1=2 



Wl 



b I- 

E 



2 — {Pr-Pj,r) 
r=l 



r=l 
2 



{Pr - Pj,'. 
PrPo 



Pr 



\w\=2 \v\=2 r=l 

We note that some of the binomial coefficients above are equal to zero. This implies that 

\2 ^ f- ^ \2 ^ I- ^ \2s 



jM\ - ^2 I, ^ PrPo ^ 



Pr 
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Using this together with 



77*1 = max —r- = max —7 > Mot. 



(see the comment after Theorem 12. ip the proof is easily completed. □ 
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